hep-ph/0608049v3 



Two Gallium data sets, spin flavour precession and 

KamLAND 

Bhag C. Chauhan*, Joao Pulido^ 

CENTRO DE FISICA TEORICA DAS PARTICULAS (CFTP) 
Departamento de Fisica, Instituto Superior Tecnico 

Av. Rovisco Pais, P-1049-001 Lisboa, Portugal 
and 

Marco Picariello-'^ 

I. N.F.N. - Lecce, and Dipartimento di Fisica, Universita di Lecce 
Via Arnesano, ex Collegio Fiorini, 1-73100 Lecce, Italia 

Abstract 

We reexamine the possibility of a time modulation of the low energy solar neu- 
trino flux which is suggested by the average decrease of the Ga data in line with our 
previous arguments. We perform two separate fits to the solar neutrino data, one 
corresponding to 'high' and the other to 'low' Ga data, associated with low and high 
solar activity respectively. We therefore consider an alternative to the conventional 
solar -|-KamLAND fitting, which allows one to explore the much wider range of the 9i2 
angle permitted by the KamLAND fitting alone. We find a solution with parameters 
Amli = 8.2 X 10~^eV'^,tan'^9 = 0.31 in which the 'high' and the 'low' Ga rates lie far 
apart and are close to their central values and is of comparable quality to the global 
best fit, where these rates lie much closer to each other. This is an indication that the 
best fit in which all solar and KamLAND data are used is not a good measure of the 
separation of the two Ga data sets, as the information from the low energy neutrino 
modulation is dissimulated in the wealth of data. Furthermore for the parameter set 
proposed one obtains an equally good fit to the KamLAND energy spectrum and an 
even better fit than the 'conventional' LMA one for the reactor antineutrino survival 
probability as measured by KamLAND. 
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1 Introduction 



After having determined that the solar neutrino problem is essentially a particle physics one 
and neutrinos oscillate [1] - [4], the next step is to search for a possible time dependence 
of the active solar neutrino flux and to investigate its low energy sector (E < 1 — 2 MeV) 
which accounts for more than 99% of the total flux. These two issues in association with 
each other may lead to further surprises in neutrino physics, possibly the hint of a sizable 
magnetic moment. Although a lot of effort has been devoted to examining the possible 
time modulation of the neutrino flux [5] - [10] , this question remains largely unsettled. The 
claim made in the early days [11, 12] of a possible anticorrelation of the Homestake event 
rate [13] with sunspot activity remained unprovcn, as no sufficient evidence was found in its 
support. More recently the Stanford Group has been claiming the existence of two peaks 
[8] in the Gallium data at 55-70 SNU and 105-115 SNU. Moreover, Gallium experiments 
[14] - [17], which have been running since 1990-91 and whose event rates are mainly due 
to pp and ^Be neutrinos (55% and 25% respectively), also show a flux decrease from their 
start until 2003 [18]. These data are hardly consistent with a constant value and exhibit a 
discrepancy of 2.4(7 between the averages of the 1991-97 and 1998-03 periods (see table I). 
No other experiment sees such variations and none is sensitive to low energy neutrinos with 
the exception of Homestake whose rate contains only 14% of ^Be. Hence this fact opens 
the possibility that low energy neutrinos may undergo a time modulation partially hidden 
in the Gallium data which may be directly connected in some non obvious way to solar 
activity [19]-[21]. Hence also the prime importance of the low energy sector investigation. 
To this end, in the near future, two experiments, Borexino [22] and KamLAND [23], will be 
monitoring the ^Be neutrinos. 



Period 


1991-97 (I) 


1998-03 (II) 


SAGE+Ga/GNO 


77.8 ±5.0 


63.3 ±3.6 


Ga/GNO only 


77.5 ± 7.7 


62.9 ±6.0 


SAGE only 


79.2 ± 8.6 


63.9 ±5.0 



Table 1 - Average rates for Ga experiments in SNU (see ref.[18]) . 

In our former work [19, 20] we developed a model where active neutrinos are partially 
converted to light sterile ones in addition to LMA conversion at times of strong magnetic 
held, thus leading to the lower Gallium event rate (II). When the fleld is weaker, LMA 
oscillations act alone and the higher rate (I) is obtained. The resonant conversion from 
active to sterile neutrinos is originated from the interaction between the magnetic moment 
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and the solar magnetic field. Its location is determined by the order of magnitude of the 
corresponding mass squared difference Am^Q and is expected to occur in the tachocline. 
where a strong time varying field is assumed. This implies Am^Q = O{10~^)eV'^ so that 
the LMA resonance and the active sterile one do not interfere ^. It was shown [19, 20] 
that the high Gallium data as in table 1 can be fitted assuming a 'quiet sun' with a weak 
field while the low data are fitted with a strong field ('active sun'). On the other hand the 
bimodal character of the data as claimed by the Stanford Group [8] cannot be explained in 
the context of this model. 

The purpose of this paper is to present oscillation and oscillation + spin fiavour preces- 
sion (SFP) fits to the data, taking Am|i, Ou as free parameters. Two separate global fits 
for Gallium sets (I) and (11) are performed with the solar data which were available in the 
corresponding periods (see table 2). 



Experiment 


Data 


Theory 


Reference 


Homestake 


2.56 ±0.16 ±0.15 


8.09±|| 


[13] 


SAGE 


see table I 


125.9±1|? 


[16] 


Gallex+GNO 


see table I 


125.9±i|? 


[15] 


Kamiokande 


2.80 ±0.19 ±0.33 


5.69 ± 1.41 


[26] 


SuperK 


2.35 ±0.02 ±0.08 


5.69 ± 1.41 


[27] 


SNO CC 


1 ao _|_0.06 _|_0.08 
i.UO 3=Q 06 =1=0 09 


5.69 ± 1.41 


[28] 


SNO ES 


9 _t_0.22 1 0.15 
Z.OO =no.22 ="=0.15 


5.69 ± 1.41 


[28] 


SNO NC 


A QA -1-0.21 _|_0.38 
4.y4 =to.2i =1=0.34 


5.69 ± 1.41 


[28] 



Table 2 - Data from the solar neutrino experiments except Ga which is given in Table 1. 
Units are SNU for Homestake and 10^cm~^s~^ for Kamiokande, SuperKamiokande and 
SNO. We use the B SOS (OP) solar standard model [24]. 

Hence for data set (1) we consider a rate fit to Gallium, Chlorine and Kamiokande 
data, the only ones existing at the time, together with the KamLAND fit, as these data are 
obviouly independent of solar activity. For set (11) we consider a global fit (rates ± spectrum) 
with the exclusion of the Chlorine one, the replacement Kamiokande SuperKamiokande 
and the inclusion of the SNO data. We consider a solar field profile peaked at the bottom 
of the convective zone. We hence determine the parameter values that lead to the best 
solar±KamLAND fit for data set (I) . Using these, we establish the best solar fit for set (II) 

^For mathematical details we refer the reader to [19, 20]. We use throughout a neutrino transition 
magnetic moment fiv = 10~^^/is- 
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on the grounds of its order parameter values Am^Q, which situates the SFP resonance, and 
Bq, the field strength at the peak. This is also the best solar +KamL AND fit for data set 
(II). The Ga (I) and (II) rate predictions lie quite close to each other in this fit, both almost 
2a away from their central values. This is not totally surprising, as their relative weight 
is small within the wealth of solar+KamLAND data. It tells us instead that the global fit 
analysis is not adequate for the investigation of the possible Ga flux variabihty. To this end 
we present a choice of a shghtly different value of Am|i , 612 whose fit to data is almost as 
good as the best fit and in which the Ga (I) and (II) rate predictions he much further apart. 
We use the BS05 (OP) solar model throughout [24]. 

Since our solar fit is independent from the conventional solar one, our only a priori 
parameter constraints come from the KamLAND fit for which the mixing angle 612 bears 
a considerable uncertainty (fig.4 of ref. [23]). In this way we scan the range tan^ ^^12 G 
[0.1,10]. Moreover in the new scenario the prediction for the KamLAND antineutrino 
survival probability is consistent with the measured one [25]. A clear distinction between 
our scenario and the conventional LMA one is expected to be provided by the forthcoming 
KamLAND data if and when new reactors start and others cease operation, thus changing 
the effective source-detector distance travelled by the antineutrinos. 



2 Global Fits and Survival Probability 

We start this section by introducing the following solar field profile 

^= cM6(f-..)] 

whose peak value Bq is situated at Xc = 0.71, x denoting the fraction of the solar radius. 
For the solar statistical analysis we use the standard function 

-1 



31,32 



where indices ji , j2 run over all solar neutrino experiments and the error matrix includes the 
cross section, astrophysical and experimental uncertainties ^. For the KamLAND analysis 
we compute the prompt energy spectrum of the positron according to 

S{E) = N j a{E,)R{E, E')^dE' . (4) 



^For mathematical details see e.g. [29]. 
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Here , E denote the physical and measured prompt event energy, so that E' includes the 
e+e" annihilation energy and the positron kinetic energy, E' — 2me + T'^+. From the well 
known relation 

Ep = mN-mp + me + = 1.804 Me^ + T^+, (5) 

(for zero neutron recoil) one gets Ep = 0.782 MeV + E' . Hence eq.(4) is also an integral 
over antineutrino energy. The quantity N is a normalization constant obtained from the 
total number of events in the absence of antineutrino disappearance above the energy cut 
at 2.6 MeV [23] 

J S{E)dE = 365.2. (6) 

The quantity a{Ep) is the total cross section for the reaction up — > ne"*" with zero neutron 
recoil energy [30] and R{E,E') is the energy resolution function 



R{E,E) = ^=exp 



{E - Ej 
2^2— 



(7) 



with a = 0.062v-E. The KamLAND energy spectrum comprises 13 bins of size 0.425 MeV 
in the range from 2.6 to 8.125 MeV. For convenience we integrate eq.(4) in antineutrino 
energy, so that in order to take into account the important low energy tail of the spectrum 
we run the integration from Ep^^^ = l.SOAMeV up to Ep^^^ = (8.125 + 0.782) MeV. In 
contrast, in eq.(6) the integral extends over all 13 energy bins. The information on the flux 
and spectra of the 20 power reactors is contained in the time averaged differential neutrino 
flux which denotes the number of neutrinos per unit energy, area and time. We have 
used the approximation (see e.g. [31]): 



d4> 



dE, 



20 n rIN 

Y.Posc{Ep,L,)ct>,Y.^^- (8) 



In this expression 



P j=l /=1 J^f d-^P 



PosciEp, L,) = 1 - sin'2esin' (^^) (9) 



is the well known oscillation survival probability formula for antineutrinos from the j*'^ 
reactor at distance Lj, (pj is the reactor flux at the KamLAND detector [32], qj and Ef are 
the relative flssion yields and flssion energies [23, 31]. Following [31], we take 



dEp 

with the coefficients given in [33]. 



exp{a} + a)Ep + ajE'^,) (10) 



4 



We have used Poisson statistics for the 13 KamLAND energy bins. Our ^ function 
therefore is [1] 

13 

Xkl — XI 

with a being an absolute normahzation constant and Osys = 6.5% the total systematic 
uncertainty. The binned spectrum calculated in the basis of eq.(4) added to the background 
spectrum [34] is thus inserted in (11) as Rf^. 

We are now in a position to obtain the result of the best fits which we show in table 3. 
For sets (I) and (II) we take fs — 1-0. We found for set (I) 

Amli = 8.2 X IQ-^eV^, tan^O = 0.31, a = 1.01 (12) 

We include for comparison the 'KamLAND only' LMA best fit [23] with all Gallium data 
replaced by their average [18] 



Raa = 68.3 ± 2.9 SNU (13) 

and with parameters ^ 





Am2i 


= 7.9 X 10-^ey^ tan^e 


= 0.46, Bo = 0, /b = 0.9, a = 


1.01. 


(14) 




Ga 


CI K (SK) SNOnc 


SNOcc SNOes xltes xk.. 


XSNO^i 


Xkl 


Set (I) 
Set (II) 


71.7 
69.6 


2.66 2.29 

2.18 5.53 


3.09 

1.54 2.16 2.28 44.6 


45.8 


15.3 
15.3 


LMA 


64.8 


2.74 2.30 5.10 


1.75 2.28 0.95 45.7 


43.1 


14.5 



Table 3 - Best fits to data sets (I), (II) and LMA best fit. For data set (I) only Ga, CI 
and Kamiokande data were available and for set (II) all SuperKamiokande and SNO data 
were available but not CI, hence the blank spaces. In set (II) only the Ca rate contributes to 
xlates- Units are SNU for Ga and CI and lQi^cm~'^s~^ for SK and SNO. 



Using Am2i, tan^O from (12), we get the best fit for set (II) with 

Am^o = -6-5 X 10-^e1/^ 2,QQkG. (15) 

^The 'LMA parameter values' we consider are the ones fixed from the best fit analysis to KamLAND 
data only, since the conventional solar fits are not taken into account in the present work. 
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r>exp 
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Tabic 3 XsKsp contains the contribution from electron scattering in SuperKamiokandc (44 
data points) and XsNOgi contains the whole SNO data set consisting of 38 data points (34 
CC, 2 NC and 2 ES day-night rates). We have 

X% = x| + Xkl (16) 

with 

Xq — Xlates + XSK^p + XSNOgf i^"^) 

For sets (I), (II) and LMA we have therefore xli/d-O-f- = 1-42, 1,15 1.12 for 13, 94 and 
93 degrees of freedom (d.o.f.) respectively. We thus see that for the best global fit the 
parameter Am^Q is well above O[10~^eK^]. This imphes that the resonances of the low 
energy neutrinos lie too deep inside the sun for the possible time variation of the field to 
significantly modulate their fiux. This is particularly true for the pp sector. In fact the 
resonances of these neutrinos lie around x = 0.45, where the field is 40% of its maximum 
and the matter density is much higher than near the peak. Hence for the best global fit 
only a small difference is expected between the two Gallium values (see table 3). 

It should be noted however that the global best fit presented here refers to 97 different 
experiments and the Gallium data account for only one of them, so their relevance is smeared 
by all other data. For a global solar fit the situation is much the same, as the number of 
experiments is reduced to 84 or 16 with the information on the low energy neutrinos included 
only in one of them. The global fit is therefore not a good measure for the investigation of 
time variability of low energy neutrinos: the information on them is nearly hidden within 
the wealth of solar and KamLAND data. 

A small change in the parameters provides a considerable difference: in table 4 we 
present the fits for Am^g = -1.7 x IQ-^eV^, fs = 1.0, fse = 1-1, a = 0.99 and Bq = 280kG 
for set (II). This leads to a much larger separation between high and low Ga rates in better 
accordance with table 1 at the price of a slightly higher x^i- From table 4 we obtain 
xli/ d.o.f. — 1.39, 1.19 for sets (I) and (II) with 13 and 94 d.o.f. respectively. All values are 
well within 2a of the experimental data. 





Ga 


CI 


K (SK) 


SNOnc 


SNOcc 


SNOes 


xLtes 


XsKsp 


XsNO 


Xkl 


Set (I) 


73.8 


2.71 


2.29 








2.80 






15.3 


Set (II) 


60.3 




2.28 


5.65 


1.59 


2.25 


0.54 


47.5 


48.5 


15.3 



Table 4 - Same as table 3 with Amfo = -1.7 x IQ-^eV^, Bq = 280kG. 



In fig.l we plot the 'KamLAND only' fit on which we superimpose our own fit (see table 
3) in the plane Am|i, tan^O. The contour lines correspond to the 95%, 99% , 99.73% CL. 
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We also show our choice of parameters (table 4) in this figure whose goodness of fit is seen 
to lie well within the 95% CL. On the other hand a comparison with fig. 4 (a) of ref. [23] 
shows that it lies inside the 95% CL of the conventional solar fit which we neglected in the 
present paper. The best fit for the 'low' Ga rate in the plane Am^Q, Bq together with the 
90%, 95% and 99% CL contours is plotted in fig. 2 where we also show the parameter choice 
of table 4 lying well within the 90% CL. 

We next describe a crucial consistency test for the present scenario, namely the survival 
probability prediction for reactor antineutrinos as measured by KamLAND [23] . Our results 
are depicted in fig.3 for the LMA+SFP fit prediction ^ with Ami^ = 8.2 x IQ-^eV^, tan^O = 
0.31 (upper curve) and for the LMA best fit, Am|i = 7.9 x lO^^eV^, tari^O = 0.46 (lower 
curve). For the average source-detector distance of 180 km, as reported by KamLAND, we 
get from fig.3, 

P = 0.576 (LMA), 0.623 (SEP fit) (18) 

to be compared with the data, P = 0.658 ± 0.064 [23]. Interestingly enough it is seen 
that the best of the two fits is the one for SEP with the parameter choice given in table 
4, namely the one leading to the two separate Gallium data sets at 73.8 and 60.3 SNU. 
The LMA solution provides the poorest fit, at 1.28(7 away from the central value. The 
present data may only allow us to expect a clear distinction between the fit of table 4 and 
the LMA one for average distances below 110-120 km (see fig.3), which will happen only if 
new reactors come into operation or others cease. Altogether the quality of the LMA and 
LMA-I-SEP energy spectra is similar {xKL/d-o.f. — 1.12 and 1.18 respectively), so that the 
two predictions are equivalent. However the accumulation of more data from KamLAND 
will no doubt clarify the situation. Besides this, we also need more data from the low energy 
solar sector to tell us whether its time modulation is a true or just an apparent effect. 

3 Summary and Concluding Remarks 

We have considered an alternative to the conventional solar neutrino data fit which takes 
into account the possible time dependence of the Gallium fiux in the form of two data sets 
corresponding to two different periods and separated by 2Aa (table 1). This is interpreted 
as an indication of a possible time variability in the low energy neutrinos, mainly pp, which, 
though constituting more than 99% of the solar ffux, only contribute with 55% of the Gallium 
flux. 

The simplest mechanism that explains the discrepancy of the two data sets is based 
on the partial conversion of active to sterile neutrinos through the spin flavour precession 

^We remind the reader that the KamLAND survival probability does not depend on Amio nor Bq. 
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originating from a varying magnetic field in addition to the LMA oscillation. In the global 
best fit that takes into account all solar data including KamLAND. the two Gallium rate 
predictions for each period lie much closer to each other than indicated by the data. Low 
energy neutrinos however contribute to only part of the Gallium rate which in turn amounts 
to one single experiment within the wealth of solar+ KamLAND data. It is therefore quite 
natural to expect that the minimum x^i is not obtained for a Gallium rate prediction close 
to its central value. For this reason the global fit analysis is not a good measure for the 
investigation of the time variability of the low energy neutrinos. We have presented a slightly 
different parameter choice with the same Arri^i and 9 values as the best LMA+SFP fit (see 
fig.l) but with Amfo = -1-7 x 10-^eV'^,Bo = 280kG (see fig.2) and for which the two 
Gallium rates are much separated from each other (see table 4). It can be seen from figs.l 
and 2 that this parameter choice lies near the edge of the 90% CL of the best fit. 

The KamLAND spectrum prediction associated with the choice of parameters of table 4 
is practically equal in quality to the LMA one as seen from a comparison between tables 3 and 
4. We also tested the antineutrino survival probability as a function of the source-detector 
distance. It was found that the parameter choice leading to the two clearly separate Gallium 
rates (table 4) provides the survival probability prediction which is in best agreement with 
the data (fig. 3). 

A close inspection of fig. 3 tells us that with the presently existing data, one can only 
hope for a distinction between the two scenarios (LMA and LMA+SFP) if the effective 
soTuxe-detector distance is reduced to less than 110 - 120 km. The Shika2 reactor, which 
started operation in late 2005 and gradually increased the operation time at the nominal 
power (4GW thermal power) from last November, will reduce this effective distance from 
160 - 190 km to 140 - 170 km [35]. Although this cannot provide a definitive answer, it 
is expected that the situation will continue to evolve as new reactors come into operation 
and others cease. The accumulation of more data from KamLAND and especially on the 
low energy solar neutrino sector is essential to disentangle the prevailing mysteries of solar 
neutrinos. 
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Figure 1: Our best fit Amli = 8.2 x 10 ^eV^, tan^O = 0.31 (table 3) is superimposed on the 
neutrino oscillation parameter allowed region from KamLAND [23]. The contour lines are 
the 95%, 99%, 99. 73% CL. 
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Figure 2: The best fit (table 3) in the plane Am^Q, Bq and the fit of table 4- The contour 
lines are the 90%, 95%, 99% CL. 
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Figure 3: The antineutrino survival probability for the fits considered: LMA (dashed line, 
table 3) and LMA+SFP ones (full line, tables 3 and 4). Note that the two LMA+SFP fits 
only differ from each other through the values of Am\Q and Bq from which the survival prob- 
ability is independent. Here the straight solid line refers to the central data point (P =0.658) 
and the two horizontal dashed lines to the la range [25]. The effective reactor distance 
(180km) is marked by the vertical dashed line. 
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